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The Telstar I satellite was launched into orbit by a three-stage Delta launch 
vehicle guided by the Bell Telephone Laboratories command guidance 
system. The Delta program is a National Aeronautics and Space Adminis- 
tration sponsored series of missile flights designed to place various scientific 
payloads into orbit around the earth. This paper discusses the theory of the 
guidance equations employed by the command guidance system in the Delta 
program. 

I. INTRODUCTION 

The Telstar I satellite was launched into orbit by a three-stage Delta 
vehicle on July 10, 1962, at the Atlantic Missile Range. The Bell Tele- 
phone Laboratories command guidance system, developed for the Air 
Force, was employed to guide the Delta vehicle. The guidance system is 
shown in Fig. 1. The missile-borne equipment, housed in the second stage 
of the Delta vehicle, serves as a radar beacon to provide return pulses 
to the tracking radar and as the receiving portion of the command data 
link between the ground and the missile. The tracking radar functions 
as the transmitting portion of the data link and also serves as a sensor 
to determine the slant range, azimuth angle, and elevation angle of the 
missile during its flight. The precision tracking radar and the missile- 
borne guidance package were manufactured by the Western Electric 
Co. The guidance computer was designed and manufactured by the 
Univac Division of the Sperry Rand Corporation. 

The three-stage Delta missile, designed by the Douglas Aircraft Com- 
pany, consists of two liquid propellant stages and a solid propellant 
third stage. The powered flight portion of the Telstar satellite trajectory 
is shown in Fig. 2. The guidance system transmits corrective pitch and 
yaw steering commands during first- and second-stage powered flight. 
Second-stage engine cutoff is ordered by the guidance system when the 
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Fig. 2 — Telstar I trajectory. 

position and velocity of the missile are such that the addition of the third 
stage velocity impulse at the end of the ballistic coast phase would yield 
the desired orbit. The unguided third stage is spin-stabilized to maintain 
attitude control. For the Telstar satellite trajectory, the third-stage 
velocity impulse was added at the perigee of the final orbit after a bal- 
listic coast phase of approximately 600 seconds between second-stage 
cutoff and third-stage ignition. The Telstar satellite was separated from 
the third stage 120 seconds after third-stage burnout. 

The guidance system steering and cutoff commands during the first 
and second stage ascent are calculated in the guidance computer, using 
the radar tracking data of the missile's position as the basic input infor- 
mation. The computer is programmed with a set of guidance equations 
that process the radar data and compute the desired commands to the 
missile. 

This paper contains a description of the theory and design of the 
guidance equations used in the Telstar satellite flight. Guidance concepts 
are presented from the point of view of orbital mechanics and control, 
followed by a description of first- and second-stage guidance. The last 
section summarizes the results achieved in the Telstar satellite flight. 



II. GUIDANCE CONCEPTS 

The Keplerian motion of an earth satellite is completely defined by 
the specification of the vector position and velocity at an epoch. The 
satellite coordinates at insertion into orbit can be expressed in terms of 
the spherical coordinate system of Fig. 3 as follows : 7 3 , 73 , and /3 3 are 
are the magnitude, the elevation angle above the local horizontal, and 
the azimuth from North, respectively, of the velocity vector; R3 , X 3 , 
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Fig. 3 — Insertion coordinates. 

and w are the radial distance from the earth's center, the longitude, 
and geocentric latitude, respectively. The elements of the ellipse, i.e., 
apogee and perigee distances, are determined by V 3 , R3 , and 73 . The 
orientation of the ellipse relative to the earth in terms of inclination, 
argument of perigee, and longitude of the ascending node depends, in 
general, on all the insertion coordinates. It would be necessary to control 
all six coordinates of the satellite as well as the time of insertion to 
achieve a specified orbit in inertial space. For earth satellites the re- 
quirements on the insertion time are usually not stringent and are largely 
determined by launch time variations. The problem of guidance thus 
consists of achieving a specified set of six insertion coordinates. 
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In the case of the Telstar satellite trajectory, the unguided third 
stage ignites at a predetermined time on the transfer ellipse following 
the completion of guidance at second-stage cutoff. Since the charac- 
teristics of the third stage are known, fixed relations exist between the 
insertion coordinates and the coordinates of the missile at second-stage 
cutoff. The transfer ellipse is denned by the position and velocity vectors, 
R2 and V 2 , respectively, of the missile at second-stage cutoff. The direc- 
tion of the velocity impulse added by the third stage is determined by 
the attitude of the missile's roll axis (axis of thrust application) at 
second stage cutoff. If we define p 2 as a unit vector lying along the roll 
axis of the missile, the eight independent coordinates of the missile at 
second-stage cutoff uniquely determine the six insertion coordinates. 
That is, 

{ V 2 , R 2 , ? 2 } -> [V z , 78, fc , Rt, Xi , ^3} • (1) 

Also the six insertion coordinates specify any six of the cutoff coordi- 
nates in terms of the remaining two. 

During the first and second stages, guidance of the missile is limited to 
steering in the pitch and yaw planes and cutting off the second-stage 
rocket engine. The missile is constrained to fly a predetermined tra- 
jectory by pitch and yaw steering during first and second stage. This 
trajectory, if followed exactly, would yield the coordinates at second- 
stage cutoff that would produce the desired orbit. However, propulsion 
system variations, deviations in the attitude control system, and radar 
noise cause dispersions in the cutoff coordinates from the expected 
values. 

As discussed in Section IV, |V 2 | and p 2 can be controlled directly at 
cutoff to provide direct control over three of the insertion coordinates. 
Control of |V 2 | by cutoff of the rocket engine provides the most sensitive 
control of V 3 . The two coordinates defining p 2 are controlled by pitch 
and yaw steering. Steering could be based on deviations of position, 
velocity, or attitude coordinates from the desired reference trajectory. 
The selection of the coordinates to be controlled by steering is deter- 
mined by the steering system design, which is based on minimizing the 
errors in the insertion coordinates. As demonstrated in Section 4.2. 
control of the pitch and yaw attitude angle, i.e., control of p 2 , affords 
the best control over insertion coordinates. 

For the Telstar satellite trajectory, |V 2 | and p 2 were used to constrain 
T" 3 , 73 , and /3 3 at insertion to provide the desired apogee altitude and 
inclination. The other orbital elements were effectively controlled by 
steering the missile to the desired reference trajectory during the first- 
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and second-stage powered flight. The equations relating |V 2 | and p 2 
to the desired orbital conditions, and the equations for pitch and yaw 
steering derived from the reference trajectory, were programmed into 
the guidance computer. The targeting task for the Telstar satellite 
mission was to determine the numerical coefficients for the equations 
used in first- and second-stage guidance. 

III. FIRST-STAGE GUIDANCE 

The missile's position as tracked by the radar in slant range, azimuth, 
and elevation angle is converted to the earth-fixed Cartesian frame 
shown in Fig. 4. The angle A is selected such that the Y-Z plane is 
approximately parallel to the pitch plane of the missile and the X axis 
lies in the yaw plane. The angle E is determined by passing the Y axis 
through the expected position of the missile at second-stage cutoff. 

From approximately 90 seconds after lift-off, pitch and yaw steering 
orders are transmitted to the missile. Yaw steering is based upon devia- 
tions in the X velocity from a reference polynomial in Y. The poly- 
nomial is selected to match the desired X component of velocity, which 
is a function of the launch azimuth, the pitch and yaw programmed 




Fig. 4 — Computational coordinate system (X axis is in the horizontal plane; 
Y axis is at azimuth angle A from true north and at elevation angle E above 
horizontal plane; Z axis is perpendicular to A' and Y axes, completing the right- 
hand Cartesian coordinate system). 
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rates in the missile, and the performance of the propulsion system. In a 
similar manner, pitch steering commands are based on deviations of Z 
from the reference velocity. The steering orders sent to the missile are 
related to the error signals in a way that balances loop response to 
missile autopilot errors, propulsion system dispersions, and radar 
tracking noise. The steering commands are transmitted to the missile 
via the radar data link at the pulse repetition rate of the radar. 

IV. SECOND-STAGE GUIDANCE 

4.1 Second-Stage Cutoff 

In this section the methods used for determining the required coor- 
dinates at cutoff will be derived. Since the Telstar satellite was inserted 
into orbit at perigee, control of apogee distance (r a ) and inclination (i) 
required that the following relations be satisfied. 



y 1 + R 3 /r a 
73 = (2b) 

/3 3 = sin -1 (cos i/cos <p%). (2c) 

As discussed in Section II, we can write 

a K =.MV 2 ,R 2 ,p 2 ), it~=lto6 (3) 

where a K represents any one of the six insertion coordinates. The six 
vector functions of (3) depend only on the transfer ellipse and the 
characteristics of the third stage. If V 2 , R 2 , and p 2 are expressed in 
terms of the coordinate system of Fig. 4, (3) can be expanded in a 
Taylor series about the expected cutoff coordinates. Linearity studies 
on the variations in the cutoff coordinates indicate that only the first- 
order terms in the expansion are significant. Equation (3) simplifies to 
the form 

^-^o = E^(C i -C 10 ) (4) 

where C,- for i = 1 to 8 represents the 8 cutoff coordinates. The partial 
derivatives in (4) are obtained by perturbing cutoff coordinates and 
integrating numerically in a digital computer through third-stage burn- 
out to determine the incremental changes in the insertion coordinates. 
Equation (2) can also be approximated by the first-order terms of a 
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Taylor expansion and combined with (4), giving three linear equations 
which can be concisely expressed by the single vector equation below. 

4,(V 2 - V 2o ) + A 2 (R - R 2o ) + A 3 ( Q2 - p 2o ) = 0. (5) 

Ai , A 2 , and A 3 are 3 X 3 matrices whose (constant) elements are 
defined by the partial derivatives used in (4) and the expansion of (2). 
Equation (5) can be solved for any three of the cutoff coordinates as a 
function of the remaining five. The three coordinates selected are 
| V 2 1 , the magnitude of the velocity vector at second-stage cutoff, and 
the unit vector p 2 as defined by the pitch and yaw Euler angles, 2 and 
fa . The missile is cut off when the measured | V 2 | satisfies (5). The 
attitude constraints on p 2 are met by comparing the values of 2 and ^ 2 
required for the solution of (5) with measured values and command- 
ing the missile to turn by the differences. 

4.2 Steering System Design 

It was shown in the previous section that if orbital elements are to be 
controlled, the relationships of (5) must be satisfied at second-stage 
cutoff. Design of the steering system is based not on minimizing the 
dispersions of the individual variables, R 2 , V 2 , and p 2 at the end of 
second stage, but rather on minimizing the orbital errors which are 
caused by errors in these variables. 

The pitch and yaw steering system design consists of finding a steering 
transfer function which minimizes insertion errors due to radar tracking 
noise and missile dispersions. Missile dispersions include propulsion 
system variations and errors in the attitude control system of the 
missile. 

4.2.1 Steering System 

A block diagram of the pitch steering system is shown in Fig. 5. A 
similar diagram could be drawn for yaw steering. 

The computer, operating on the radar data, obtains its measure of 
the vehicle position in the Z direction, Z c . After steering has started, 
the measured trajectory variables are compared with a reference tra- 
jectory and corrective pitch turning rates, d„ , are sent to the missile. 
The pitch rate programmed in the missile, P , and the ordered turning 
rates are the inputs to the autopilot's reference integrating gyro. The 
function of the autopilot is to align the direction of the acceleration 
vector, 0, with the desired attitude, as indicated by the gyro output, 6 T . 

If the missile's roll axis is aligned with the Y axis of Fig. 4, a small 
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Fig. 5 — Pitch steering system (a is the missile thrust acceleration; 6 P is the 
programmed turning rate; Eg p is the error in the programmed rate; Eg is the error 
in the thrust vector direction; Ey. is the radar error in measuring the position Z; 
Z c is the computer's measurement of Z; and B isthe guidance ordered turning rate). 



change in pitch attitude, 0, multiplied by the thrust acceleration, a, is 
the change in Z-direction acceleration, Z. Two integrations, represented 
by their Laplace notation in Fig. 5, then give the vehicle position, Z. 

Some of the missile error sources are gyro input errors, propulsion 
system errors, and misalignments between the direction of the vehicle 
acceleration, 8, and the reference attitude, 6 T . The gyro input errors are 
caused by errors in the programmed rate and by electrical or mechanical 
unbalances. Because the reference trajectory is determined by the pro- 
grammed turning rates and the expected performance of the propulsion 
system, propulsion system variations can be replaced in the block dia- 
gram by equivalent programmed rate errors. All of the attitude errors, 
including the integrated rate errors, are caused by slowly varying dis- 
turbances. 

The radar noise error in Z, E z , is the product of the elevation angle 
error and the distance from the radar to the missile. The radar noise 
power spectrum has most of its energy at frequencies higher than those 
encountered in the missile attitude disturbances — a fact important in 
the steering loop optimization. 
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Fig. 6 — Simplified steering loop 
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Fig. 7 — Velocity errors. 

4.2.2 Steering Loop Optimization 

In Fig. 6, the programmed pitch rate and the reference trajectory it 
describes have been removed from the steering loop. All of the missile 
dispersions are combined in a single attitude error source, E 6 . With the 
reference trajectory removed, 0, Z, Z, and Z represent dispersions about 
the reference values. The quantity to be minimized by steering may be a 
linear combination of some of the above variables, as will now be shown. 

Fig. 7 shows the effect of 6 and Z on the total insertion velocity vector 
of the Delta missile with its unguided third stage (assuming 6 is small) 
to be 



= Z+ V t 



(0) 



The subscript in Fig. 7 indicates reference values, and V 1U is the mag- 
nitude of the velocity increment of the third stage. In complex frequency 
notation 



M = - 



S \ W2/ 



-UW = 0(8) Y.,e 



(7) 



where o> 2 = a/V JU , and in general Y a/b is the transfer function from b to 
a. The total error in terms of the attitude and noise errors can be ex- 
pressed as 



e(s) = Y t/Ez E z +(Y t/e + ^f^JEe. 



(8) 
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The total variance of the error, a t 2 , is the integral over all frequencies 
of the power spectrum of the error signal. The optimization problem is 
to determine the transfer function, Y (lEz , (which in turn determines the 
steering equation, Y,, , Fig. 6) that minimizes a ( 2 . The fact that the 
spectral density of the attitude errors is concentrated at very low fre- 
quencies while the radar error power contains higher frequency com- 
ponents suggests that the error may be minimized by a frequency- 
selective steering equation. 

The attitude errors can be approximated by a Markov power spec- 
trum, having a low cutoff frequency, 

P„(a>) - ^ali" (9) 

where coo is very small. Using techniques described by Bode and Shannon 1 
for minimizing the mean squared error, the optimum steering equation 



is: 



YL ~ 2T* + 2 -T + (t> + ^ ) ' (l0! 



s -\ s 

W2 



where 



(2 \ 1/6 

5T-T ) • (ID 

a- o-gfl co8 co K / 

For the Delta second stage, l/o> 2 equals F m /a. If typical numbers 
for Fin , a, and T are inserted, then l/cu 2 » T, which leads to the approxi- 
mation : 



Yi„ = I T 2 \ (12) 

''II cyry I rp„ i ' \ ' - 



r (i + ts + 1 2 S A • 



The two differentiations of the position data indicated in (12), together 
with the division by a, Fig. 6, convert the position data to attitude data 
which is in turn smoothed to give the ordered turning rates. 

In the guidance equations the pitch and yaw attitude of the missile 
is determined by operating on the position data of the missile. The 
total acceleration of the missile is computed by taking second differences 
of the position data expressed in the coordinate system of Fig. 4. 
Gravity, Coriolis and centripetal accelerations are subtracted from the 
total acceleration to obtain the thrust component of acceleration. Since 
the thrust acceleration vector is aligned with the missile's roll axis and 
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the missile is roll stabilized, the pitch and yaw attitude angles can be 
computed. These are compared with the desired values to obtain pitch 
and yaw attitude errors. The desired attitude angles are computed 
from polynomials in time designed to match the programmed rates 
built into the missile's attitude control system. 

In this paper, all transfer functions have been given in continuous 
form (as functions of s), although the transfer functions are actually 
converted to digital form (functions of 2 = c' T ) before being programmed 
into the digital computer. This conversion does not significantly change 
any of the relations given here, because the frequency ranges important 
in the above equations are much less than half of the sampling frequency. 

4.2.3 Accuracy Improvement by Final Value Control 

Equation (12) is the optimum steering equation for minimizing the 
mean squared error. However, the injection errors are functions of the 
error at one instant of time, second-stage cutoff. Examination of the 
system's response suggests means of improving the steering design. 

The slowly varying attitude disturbances represented by a Markov 
distribution may also be represented by initial attitude and velocity 
errors and a constant gyro drift. The system's response to these errors, 
derived from (12), is shown in Fig. 8 as a function of time normalized 
with respect to the steering parameter T. 

Since the Delta design is not restricted to a continuous control system, 
velocity errors at second-stage cutoff in excess of those obtained with 
the optimum continuous design can be allowed. These velocity dis- 
persions can be measured and used to aim the third stage in accordance 
with (5). The velocity errors shown in Fig. 8 are either constant or 
linearly increasing with time and can be measured by filtering the posi- 
tion data. The steady-state attitude error response to gyro drift is a 
constant. If gyro drift is a major attitude error source, continuous closed 
loop control of attitude may also be relaxed, and the constant attitude 
error can be measured with a polynomial filter. 2 The errors can be 
corrected just prior to cutoff by turning the missile. The reason for 
relaxing the attitude and velocity dynamic control [making T greater 
than the "optimum" value in (11)] is that doing so decreases the atti- 
tude error caused by the effect of radar noise on the steering system. 

The outputs of the attitude and velocity filters used to measure the 
dynamic residuals also have errors caused by the tracking noise. These 
decrease as the filter length, T F , is increased. Since the attitude error 
signals involve second derivatives of tracking data, the attitude filters 
require long smoothing times. However, if the filters are started too 
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Fig. 8 — Dynamic steering response: (a) response to initial attitude error; (b) 
response to initial velocity error; and (c) response to gyro drift. 



early in the second stage, they are subject to increased dynamic error 
because the steering system has not yet settled to the steady state 
(Fig. 8). Thus, optimum filter design is dependent on the steering 
parameter T, 

For the Delta missile, the critical design criterion was that of attitude 
control. That is, for a system with optimum attitude control, the opti- 
mum velocity filter has negligible error, and the measured dispersions 
can be used as a basis for an attitude correction (aiming the third stage). 

For given dynamic error sources and filter length, the total attitude 
error, after correction of the measured error, is the combination of the 
filter's dynamic error and the noise error of the filter plus the steering 
noise error, Fig. 9. The steering system's noise error is a function of T 
and decreases as T increases. The filter's noise error is independent of 
T but decreases as the filter length, T F , increases. The filter's dynamic 
error is a function of both T F and the transient response of the steering 
system as determined by T. Because the steering and filter noise errors 
are highly correlated, they are added together directly and root sum 
squared with the filter's dynamic error. T and T F are chosen from Fig. 9 
to minimize the total attitude error. 
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Fig. 9 — Final attitude errors as functions of T and T/ . 

In summary, the Delta second-stage guidance equations give closed- 
loop control of missile attitude. Relatively large dynamic errors in 
velocity and attitude are permitted in order to decrease noise errors, 
because the dynamic errors can be measured and corrected. An attitude 
order to correct the measured attitude error and compensate for velocity 
and position dispersions in accordance with (5) is sent to the missile 
shortly before second-stage cutoff. 



V. TELSTAR I SATELLITE FLIGHT RESULTS 

The steering history of the first and second stage of the Telstar I 
satellite flight is given in Figs. 10 and 11. The maximum steering orders 
during first-stage guidance were 3.8 degrees in pitch and 1.6 degrees in 
yaw. Steering in second-stage guidance reached a peak of 0.25 degree in 
pitch and 1 degree in yaw. Corrective commands of 0.35 and 0.13 degree 
were issued in pitch and yaw, respectively, just prior to second-stage 
cutoff. The second-stage engine was cut off at 269.6 seconds from liftoff, 
as compared to a preflight value based on average propulsion perform- 
ance of 273.3 seconds. The orbital results as determined by NASA are 
given in Table I with the preflight reference values. 
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Fig. 10 — First-stage steering. 
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Table I — Orbital Results 





Pretlight Reference 


NASA Minitrack 


Altitude 






Apogee (nm) 


3000 


3044 


Perigee (nm) 


500 


515 


Eccentricity 


0.2407 


0.242 


Period (minutes) 


156.47 


157.82 


Inclination (degrees) 


45 


44.79 


Argument of perigee (degrees) 


166.6 


165 


Right ascension of ascending 


203.0 


204 


node (degrees) 
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